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Motivated by the computation of scattering amplitudes at strong coupling, we consider 
minimal area surfaces in AdS^ which end on a null polygonal contour at the boundary. We 
map the classical problem of finding the surface into an SU{4) Hitchin system. The polygon 
with six edges is the first non-trivial example. For this case, we write an integral equation 
which determines the function of the shape of the polygon. The equations are 

identical to those of the Thermodynamics Bethe Ansatz. Moreover, the area is given by 
the free energy of this TBA system. The high temperature limit of the TEA system can be 
exactly solved. It leads to an explicit expression for a special class of hexagonal contours. 



1. Introduction 



In this paper we consider a geometrical problem involving minimal surfaces in AdS. 
We prescribe a polygonal contour on the AdS boundary and we consider a minimal surface 
which ends on this contour. We devise a method for computing the area of the surface as 
a function of the boundary contour. 

It is a generalization of the problem of finding the shape of soap bubbles, or Plateau 
problem to AdS space. In the 19th century the study of surface tension was interesting 
for what it could reveal about interfaces and the size of atoms. For us, the study of this AdS 
"soap bubbles" is interesting for what it could eventually teach us about their constituents, 
gluons, and their scattering for all values of the coupling. While this is our motivation, the 
computations in this paper will be restricted to classical geometry (or strong coupling). 

1.1. Gauge theory motivation 

Our main motivation for focusing on this problem comes from the study of Wilson 
loops or scattering amplitudes in = 4 super Yang Mills theory. These observables are 
interesting in general gauge theories and contain a wealth of dynamical information. Since 
planar Af = 4 super Yang Mills theory is integrable 0,^ , one hopes to be able to compute 
them at all values of the coupling. At strong coupling they can be computed in terms of 
classical strings in AdS^ [Q. The classical equations of motion for a string say that the 
area should be extremized. Thus, we end up studying minimal surfaces in AdS^ space. 
This problem becomes tractable because it is classically integrable P,^. In fact, we will 
use its classical integrability to find equations that determine the area. 

We hope that our classical analysis will be a useful starting point for solving the 
problem at all values of the coupling. In fact, the final classical problem has a structure 
that looks like the Y system for the Thermodynamic Bethe Ansatz. 

For the computation of operator dimensions the solutions to the classical problem 
1^,0 were very useful for eventually determining the full quantum solution [^^,|T^,|TT| . We 
hope that the same will be true for the computation of scattering amplitudes. 
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Fig. 1: We have a polygonal loop with null sides on the boundary of AdS. This 
polygon is denoted here by the jagged lines. We consider a minimal surface in the 
bulk that is ending on the polygon. 



We consider surfaces that end at the AdS boundary on a special contour which is made 
out of nuU segments, see fig. 1. We are given a sequence of n points on the boundary, such 
that consecutive points are at a nuU distance, xf ^j^^ = 0. This defines a polygon with null 
sides. Conformal transformations act on this polygon and can change the positions of the 
points. One can define conformal invariant cross ratios which do not change under con- 
formal transformations. These characterize the polygons up to conformal transformations. 
We will choose the points in such a way that the distance between all points which are 
not consecutive is spacelike. In this case, the minimal surface that ends on the contour is 
spacelike. It is a spacelike surface embedded in the Lorentzian signature AdS^ space. Its 
area is infinite. This infinity is due to the IR divergencies of amplitudes or the UV diver- 
gencies of Wilson loops. One can define a renormalized area which is finite and conformal 
invariant. Our goal is to compute the renormalized area as a function of the cross ratios. 
In the context of scattering amplitudes, this quantity is called "remainder function" . See 
1^,0 for a perturbative, two loop evaluation of this function for six particles. 



1.2. Outline of the method 

The classical action for a string in AdS is conformal invariant. Thus, it contains a con- 
served holomorphic stress tensor T{z), T{z) on the worldhseet. The Virasoro constraints 
are setting these to zero. In this case, we also have spin four holomorphic quantities P{z), 
P{z). These are holomorphic and non- vanishing after we impose the equations of motion 
and the Virasoro constraints. These are additional local conserved quantities associated 
to the integrability of the theory. These currents play a crucial role in our method for 
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solving the problem. Via worldsheet conformal transformations z ^ w = f{z) we can 
locally set them to one P{z) — > P{w) = 1. This however, cannot be done globally without 
introducing some singularities. In our case the worlsheet can be parametrized by the whole 
complex plane and P{z) is a polynomial. The degree of the polynomial determines the 
number of sides of the polygon. The coefficients of the polynomial (together with some 
extra parameters we will discuss later) encode the shape of the polygon . 

The problem is integrable and it has a one parameter family of fiat connections 
parametrized by a spectral parameter ^. This family is encoding the infinite number 
of conserved non-local charges of the theory. 

The flat connection goes to a constant at infinity and the solutions of the linear 
problem {d + A)ip = grow at infinity. The solutions of the linear problem, for C = 1, 
encode the AdS coordinates and the fact that they go to infinity as 2; — > oo implies that 
they go to the boundary of AdS. They go to different points on the boundary of AdS in 
different angular sectors of the z plane. These sectors are called Stokes sectors. In any 
problem which has Stokes sectors we can define cross ratios, which end up being the same 
as the spacetime cross ratios when C = 1- 

The advantage of introducing the spectral parameter 7^ 1 is that the problem of 
determining the cross ratios as a function of the polynomial coefficients simplifies when 
C ^ or when ^ — )■ 00. In this limit we can determine the form of the cross ratios as a 
function of the coefficients of the polynomial P{z). In addition, we expect that the cross 
ratios are analytic functions of ( away from ( = 0,oo. The cross ratios as a function of ( 
display a sort of Stokes phenomenon as C — > 0, 00. One can define certain cross ratios which 
are analytic within certain angular sectors in the plane and have simple asymptotics for 
C — )• 0,00. When we change angular sectors these cross ratios have discontinuities which 
are simple functions of the cross ratios themselves. Patching together the different angular 
sectors one defines a Riemann Hilbert problem which consists in finding the corresponding 
analytic functions in each sector with the specified jumps between sectors. This problem 
can also be rewritten as an integral equation for the cross ratios as a function of C. 

Intriguingly, this integral equation looks very much like a Thermodynamic Bethe 
Ansatz equation (TBA), or Y-system, for some particles whose physical significance is not 
clear to us. Moreover, the renormalized area can be expressed in terms of the cross ratios 
via a formula which is the same as the formula for the free energy of the Thermodynamic 
Bethe Ansatz. 
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This paper is organized as follows. In section two we recall some general properties 
of strings in AdS^. We review the Pohlmeyer reduction and the associated Hitchin system 
and flat connection. We explain how the linear problem for the flat connection lets us 
recover the surface. In section three we explain how to derive the TBA equation from 
the Hitchin system. First we use a WKB approximation to understand the behavior of 
the cross ratios as a function of ( near C = 0, oo. With this information, we write down 
an integral equation. Finally we discuss some general aspects of the kinematics of the 
six-sided polygon. 

The connection between classical differential equations and the TBA system for con- 
formal quantum integrable theories was studied before in [1I^JT5| , |T6| . In fact, our results 
can be viewed as the extension of those results for massive integrable theories. 

In section four we show that the area can be computed from the free energy of the 
TBA equations. We also discuss a family of exact solutions of the TBA system which 
corresponds to the large temperature limit, in the TBA language. 

Finally, we present our conclusions. 



2. Classical Strings in AdS^ 

2.1. Pohlmeyer reduction and integrability 

Classical strings in AdS spaces can be described by a reduced model which depends 
only on physical degrees of freedom [1T7| , |T8| , |T9i^ , [2l| , ^ |23|| . We use embedding coordinates 



Y with = —1. The conformal gauge equations of motion and the Virasoro constraints 
are 

ddY - {dY.dY)Y = , dY.dY = BY.BY = (2.1) 

where we parametrize the (spacelike) world-sheet in terms of complex variables z and z. We 
can form an orthonormal basis q = {qq, ■ ■ ■ , 1J5) with qo = y, 54 = e~"/^9y, ^5 = e~°'^'^dY 
and qi = Bj where / runs from 1 to 3. 

= dY.dY, Bj.Y = Bj.dY = Bj.dY = 0, Bj.Bj = hu (2.2) 

The other inner products are flxed by the Virasoro constraints and the fact that Y.Y = —1. 
The real vectors S/'s span the normal space to the world-sheet inside AdS^. We pick 
them to be orthonormal. The signature of the inner product hjj is (—1,1,1) for standard 
Wilson loops in the usual (3, 1) signature, (—1, —1, 1) for Wilson loops in (2, 2) signature 
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and (—1, —1, —1) for Wilson loops in (1,3) signature0. There is a residual gauge freedom, 



either SO{l,2), SO{2,l) or SO{3) that rotate the B's, which we don't fix, see |2^] for 
gauge fixed equations. 

i?3 can be set to a constant, reducing the problem to AdS4. Setting both S3 and B2 
to a constant reduces the problem to AdS^. Given these vectors we can further define 



vi = Bi.d^Y, vi = Bi.d^Y 



(2.3) 



Note that with our conventions, as the B's are real, v and v are complex conjugate. We 
can define an SO{2, 4) (or SO{3, 3), or SO{4, 2), depending on the signature of spacetime) 
fiat connection on the worldsheet describing the parallel transport of q. 



dq = -Azq, 



dq = -Asq 



(2.4) 



Explicitly, 
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(2.5) 



Here djj = dBj.Bj and we lower and raise indices with h. Az is given by a similar 
expression obtained by interchanging v ^ v, djj — )■ djj, permuting the last two rows and 
permuting the last two columns. The fiatness condition for the connection A is equivalent 
to the following equations 



(2.6) 



dda - - e'^v^vi = 0, 
dvi — d/vj = 0, Bvi — d/vj = 0, 
{viv-^ -viv^) = f/, f = dd-dd+[d,d] 



f is the field strength of the gauge field d. These follow from the original equations ( |2.1| ). 
Conversely, one can first solve these equations, (|2.6|), and determine the connection A. We 



^ These changes of signature correspond to the standard procedure of analytic continuation 
in the external momenta. In the field theory language, we are changing the external momenta 
but we are not changing the usual Feynman prescription for the internal momenta, in contrast to 



what was done in |24]. In particular, the IR divergencies have the same form in all signatures. 
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can then find six independent solutions of tlie linear problem (2]1). We can orthonormalize 
these solutions at some point. Then these six solutions will remain orthornomal throughout 
the worldsheet. We have now an orthogonal matrix qf^ where A labels the solution number 
and it is the spacetime index that all vectors qi carry. The first row of the matrix, 
Qo = Qo ~ solves the original equations of motion and Virasoro constraints. The area 
is then given by 

A^2 [ d^ze" (2.7) 



This is the area in units where the radius of AdS is set to onea. Then the amplitude or 
Wilson loop is given by 

Amplitude ~ (W) ~ e"i^^ ; - — - = (2.8) 

where we also have given the expression for the radius in terms of the 't Hooft coupling 
for = 4 SYM. For other theories we have other relations between the radius and the 
microscopic parameters, but our result will continue to be the leading approximation when 
R^/a' > 1. 

An immediate consequence of (|2.6|) is that 



d{v^vi) = ^ v^vi = P{z) (2.9) 

with P{z) some holomorphic function. Since g is a complete basis, we can write d^Y in 
terms its elements 

d'^Y = dadY + v^Bj P{z) = d^Y.d^Y, P{z) = d^Y.d^Y (2.10) 

which gives a simple expression for this holomorphic function, and his anti-holomorphic 
companion, in terms of space-time quantities. Actually, one can verify directly from ( |2.1J ) 
that P defined as P = d^Y.d^Y is holomorphic. We can locally define a new variable w 
through 

dw = P{zf'^dz (2.11) 

This can be viewed as a worldsheet conformal transformation z w{z). In the new 
coordinates P{w) = (Jf;)^-P(^) = 1- This change of coordinates is locally well defined. 
However, at points where P has a zero, we have a branch cut in the w plane. We can view 



We have that j d^z = J dxdy, with z = x + iy. 
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the w plane as a sophisticated hght-cone gauge choice. This w "plane" is useful to think 
about the asymptotic structure of the solution. 

The integrability of the problem becomes manifest if we decompose A in two parts, 
A = A + where A rotates {Y, Bj) and {dY, dY) separately among themselves, and $ 
mixes them. 
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(2.12) 
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(2.13) 





where the missing elements of the matrices are all zero. It is easy to show that D = d+ [A 
and ^» satisfy the Hitchin equations: 

D,^, = 0, = 0, 

A simple consequence is that the spectral connection 

Wi = D, + C'^z, Vi = D,- + C$.- (2.14) 

is flat for all values of the spectral parameter Notice that unless C is a phase, the 
connection lies in SO {6, C). 

One can solve the linear problem V^g[C] = to define a whole family of solutions q[Q. 
It is easy to verify that Y[(] — qo[(] is a complex solution of the usual equations of motion 
and Virasoro constraints. If is a phase, Y[i^] can taken to be real. 

It turns out to be useful to study the action of the connection A on spinors of 5*0(2, 4) 
(or 5*0(3, 3), or 5*0(4, 2)). In practice, we can pick a set of gamma matrices and use them 
to rewrite ^ as a 4 x 4 matrix. With slight abuse of notation we will use the same symbol 
for the 6x6 and the 4x4 connections. We can keep the decomposition of A into A and 
$ simple, by using a set of gamma matrices adapted to their block-diagonal form: 

-l/2a J \ 

(2.15) 







-i'3 



V2 



-VjT 



V2 







1 / -da + dijT^J 
4 V da + dijT^J 



Here are appropriate Pauli matrices for 5'0(1,2), S'0(2, 1) or 5'0(3). Furthermore, 
$2 = and Az = —Az, where the overhne indicates hermitian conjugation, together 
with conjugation by either or r^^ for the SO{l, 2) or SO{2, 1) cases. 

Although $ and A satisfy Hitchin's equations, they are not the most general solution. 
There is a simple constraint which can be added to Hitchin's equations, such that all 
solutions satisfying the constraint will be of the form ( |2.12| ). If we introduce the constant 
matrix 

C=i!i (2.16) 



za2 

then the solutions of CA^C~^ = -A and C$'^C~^ = z$ have the form (|2.15| ). Hence 
classical string solutions in AdS^ correspond to solutions of Hitchin's equations which are 
fixed by the Z4 automorphism A -CA^C'^ and $ -zC$'^C~^. 

In terms of the vector- valued connections, we can construct a 6 x 6 matrix C combining 
a reflection of qi, i = 1, 2, 3, and a rotation of 7r/2, — )■ iq4, q^ — )■ —iq^, so that the Z4 
automorphism A — )■ CAC~^ and $ — > zC$C~^ has fixed points of the form ( |2.12| ) only. 
In particular, the symmetry gives V*^ = CV'C"^, and Y[Q must be related to Y[iQ] by a 
space-time rotation. 

2.2. Boundary conditions for the connection at infinity 

In this subsection we discuss the boundary conditions at infinity for the different 



components of the fiat connection. We start with the simpler case of AdS^ and then 
express the boundary conditions in such a way that they are easy to generalize for the 
most general case. 

The reduction to AdS2, amounts to setting B2 and S3 to a constant. Hence = 
0, f 2 = f3 = and vl = —P{z) vi = iP^I"^ . It turns out that in the AdSj, case, P is the 
square of another polynomial P oc . Writing the four by four connections in terms 

of two by two blocks we obtain 



A. = ^dai . 1, = 1 ^ I, A, = -Al = $+ 

(2.17) 



Flatness condition implies the following equation for a 

dda-e'' + e-'^\P{z)\=0 (2.18) 
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This generalized sinh-Gordon equation can be brought to a more standard form by defining 
a new variable w such that dw = P{z)^^^dz, (|2.11|) , and making a field redefinition 



a = a+^\ogP{z)P{z) (2.19) 

As argued in [|6[^, the appropriate boundary conditions are that a vanishes at infinity, 
decaying exponentially. One simple way to understand this is to notice that the four cusp 
solution is simply d = and P = 1. After we go to the w plane, this is indeed the 
asymptotic behavior at large w for any polynomial P. This is what we expect, since near 
each cusp the solution should be very similar to the four cusp solution. In the four cusp 
solution the w plane and the z plane coincide and we have four cusps corresponding to the 
four quadrants of the w plane. For a polynomial of degree which goes as P ~ + . . . 

Ar+4 

we have that w ~ z^^ . This implies that going once around the z plane we go around 
the w plane ^^^^ times. This implies that the total number of quadrants in the w plane 
is n = A^ + 4. This is the total number of cusps. In conclusion, a polynomial of degree A^ 
leads to an n = A^ + 4 sided polygon. 

Such boundary conditions can be understood in the following alternative way. If we 
set a = 0, then we can diagonalize $ and A by an honest gauge transformation h (the 
exponential of an anti-hermitian matrix) 



-iPizY/^ 

-p{zY/^ 

\ iP{zY'^J 
h~^A^h + h~^dh = 



(2.20) 



The statement that d — )■ at infinity is equivalent to saying that exists a gauge transfor- 
mation such that 

at infinity. We have similar expressions for with the same h. 

The boundary conditions at infinity for the general case of AdS^ can be stated in a 
similar way. After all, near each cusp we expect the same behavior as the one for the four 
cusp solution. In fact, these are standard boundary conditions from the point of view of 
the Hitchin system. There exist a gauge transformation such that h~^^zh ^diag and 
h~^^zh — )■ ^*diag infinity. Notice that the boundary conditions h~^^zh — )■ ^*diag involve 
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a simple complex conjugation. We also have that d — )■ at infinity, where a is given by 



For the case of AdS-^ and AdS^ the matrix dijT^"^ is diagonal, and has an off-diagonal 
component only for the case of AdS^. Next, we consider A^ in the same basis that diag- 
onalizes $2 at infinity. Flatness of $2 requires A^ to be diagonal as well. There are two 
further requirements: in the original basis A^ should be block-diagonal, and single- valued. 
This imposes the following boundary conditions 



where m is minus the conjugate of m for (3, 1) or (1,3) signature and plus the conjugate 
of m for (2, 2) signature. The constant m is allowed to be non-zero only for a Pi^z) of 
even degree n: in the diagonal basis, ^diag comes back to itself under z — )■ e^'^'^z up to 
conjugation by a shift matrix. In that basis, our asymptotic choice for A is invariant under 
that shift only for even n. Note that the Z4 action that restricts the form of the general 
Hitchin equation to our form (discussed around ( |2.16D ) has a simple action in the basis 
that diagonalizes $ ( |2.20| ) . The matrix C simply performs a cyclic permutation of all four 
eigenvalues. It also permutes the eigenvalues of A. Thus the condition A = —CA^C~^ 
implies that A has the matrix structure in (|2.21|) . 

As already mentioned, for the cases of AdSs and AdS^, the gauge field is diagonal 
before the gauge transformation. In particular, its structure forbids the extra term in 
( p.21|) . This term is present only for the case of AdS^. This is related to the fact that, 
for AdSs^4, once we choose a polynomial P{z), there is a unique solution satisfying the 
appropriate boundary conditions, but this is not the case for AdS^,. 

We can use these observations also to count the number of parameters that we expect 
in general. In order to count the moduli of the Hitchin equations we can go to the region 
where the zeros of the polynomial are widely spaced. In this region the matrix $ is diagonal 
in most of the space. Let us start by counting the number of parameters in the polynomial 
P. We can always set the coefficient of the highest power to one by a rescaling of z. We 
can also set the coefficient of the next term to zero by a translation. Thus we have only 
— 1 complex coefficients left. As we explained above, as long as we can set d ~ we 
can also diagonalize A. We will then have various Wilson lines for the connection around 
various one cycles of the surface. In order to have a non-zero value for the Wilson line 
we need a one cycle that is consistent with the projection. Instead of counting one 



(ill)- 




(2.21) 
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cycles we can count differential forms which are odd under the generator of Z4. We have 
a curve x'^ = P{z) with the differential xdz. The Z4 projection on the Hitchin system 
maps X — )■ IX. The one forms can be written as ry = q{z)^. For large z we require that 
they decay as j2dz or faster. The one form decaying as I/2 was already included in the 
discussion around ( |2.21D . For a polynomial of even degree we have q ~ ;^n/2-2 _(_... 
We have N/2 — 1 complex coefficients!. The total number of real parameters is then 
2(A^ - 1) + 1 + 2(A^/2 - 1) = 3(A^ - 1) = 3n - 15 which is the expected number of 
independent cross ratios for an n sided polygon (recall that n = N + 4). 

2.3. Behavior at infinity of the solutions of the linear problem 

Here we show that the above conditions on the connection imply that the solutions of 
the linear problem are such that the surface ends on the boundary of AdS on a polygonal 
contour with light-like sides. 

The above boundary conditions imply that the behavior at infinity of the fiat section 
q is controlled by the degree of the polynomial P. The large \z\ region is divided into 
71 = + 4 angular sectors of equal width. In the sector Vi the six components of Y grow 
exponentially with a specific fixed ratio, so that the worldsheet reaches a specific point Xi 
at the boundary, where the i-th cusp resides. Alternatively, Y ~ yiexpS'(z) where S{z) 
grows like a power of |z| and yi is a null vector, which identifies a point on the boundary 
of AdS^. 

Here we can think of the boundary of AdS^ as the set of null vectors y G B?'^ with 
y"^ = 0, together with the identification y ~ \y. Thus, only the direction of the null vector 
y is important^. 

As we change sectors, the point Xi (or yi) on the boundary jumps, in a way which 
is well described by Stokes theory for the fiat connection A. It is useful to consider the 
spinor version of the connection. Let ipa be a set of four linearly independent spinor fiat 
sections satisfying 

dtl^a = -A.i^at dtl^a = 'A^i^a (2.22) 



^ The fact that the dz/z term counts as only one real coefficient is related to the fact that we 
just have one compact cycle associated to this one form, which is a circle around infinity. 

^ More explicitly, we can say y = {y+,y-,yp,), /U = 0, f,2,3. = y+y_ +2/^2//^ = 0- Then 
= ^ are the ordinary Poincare coordinates. Also, given two points y and y we have that 
2y.y = —yjf-yjf-{x — xY . Thus, the inner product of two y vectors is proportional to the distance in 
R^'^ . If we form cross ratios the proportionality factors cancel. Thus cross ratios can be written 
in terms of ratios of products of projective y vectors. 
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We can write a general solution to the vector-valued linear problem as 



qflAah = tpaa^'^^^ph 



(2.23) 



Here are the six space-time gamma matrices, a, h are space-time spinor indices. Vj are 
the six gamma matrices used to convert the connection A from the vector to the spinor 
representation, a, (3 are internal spinor indices. These are the indices acted upon by the 
connection. Both gamma matrices are antisymmetric in their indices. 

The boundary conditions described above imply the existence of n sectors Wi (dis- 
placed by ^ with respect to the Vi) in which the ip"" grow exponentially at large 1^1 . Much 
like in the vector case, this determines a direction (a "twistor") Xa^i in each sector, as 
V'aa ~ Aa^jCa^i exp ^^(z), where Cq, is the direction of this growing solution in the internal 
space. The location of the space-time cusp associated to a sector Vi is determined by the 
twistors for the sectors Wi and VFi+i overlapping with Vi 



all 



(2.24) 



To appears here because the component go = is the solution Y^. The spinors \a,i are 
determining the direction of the vector yf". The factor in parenthesis is simply an overall 
scale. Thus, we can assign a spinor to each side of the polygon. The vertices of the 
polygon are then given by Ui ~ AiAi_(_i. These spinors are the same as the momentum 
twistors introduced in . 





(a) 



(b) 



Fig. 2: In (a) we see the polygon labeled by the position of the vertices Xi, or 
equivalently, the momenta pi. In (b) we see how the data defining the polygon 
emerges in our case. We have a spinor of SU{2\2) ip associated to each edge. 
The positions of the vertices yi are bilinears of spinors. These coincide with the 
momentum twistors introduced in 
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The inner product can be written as I 



yy' = {y laM ■ 7cd)e"'^' = det(AiA2A;Ay (2.25) 

where Ai^2 are the two spinors producing y and A'^^ 2 are the two spinors producing y' . It 
is then clear that not only yf = 0, but also yi ■ yi^i = = {xi — Xi^i)"^, hence the sides of 
the boundary polygon are null, as desired. 

2.4- Linear problem at infinity and Stokes phenomena 

In this subsection we describe the nature of the Stokes phenomenon at large z. We 
also express the Stokes data in terms of the spacetime cross ratios. 

Once we have understood the asymptotic form of the connection, it is straightforward 
to understand the asymptotic form of the solutions t/^a of the linear problem ( |2.22| ) . As 
already mentioned, we can perform a gauge transformation in such a way that both $2 
and $^ become diagonal at infinity. In such basis, a generic solution takes the asymptotic 
form 

i; ^ eiz"'z'^e^^'"+''^ +e2Z-^z-'^e-^^'"-''^ +esz^z'^e-^^''+''^ +e^z-'^z^^ 

(2.26) 

where w = J P{zY^'^dz ^ z^ , w = J P{z)^^'^dz ~ z^ . This asymptotic solution is a rather 
reliable approximation of the exact solution at large \z\, as long as one restricts the range 
of arg z appropriately. The failure of the approximation is governed by Stokes phenomena. 
Consider the behavior of a generic approximate solution as one varies argz at large \z\: 
the four exponential terms will take turns at controlling the asymptotic behavior, so that 
t/j will be proportional to, say, ei whenever Rew is positive and bigger than \Imw\ so that 
{w + w) is the largest exponent. A generic exact solution will have a similar behavior, and 
in each sector 

2n 3n 2n n 

Wi : — i < arg z < — i 

n n n n 

it will point in some direction Ai, but At+4 does not have to coincide with A^. We will 
denote as 5"^ the appropriate largest exponent. 

^ The determinant of (Ai,A2, A3A4) is equal to e"'^'''*AQ,4A/3,2A-y,3A5,4- 
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Fig. 3: The light hnes show the four exponents of the WKB approximation for 

large I X I clS function of the argument of z. The darker line follows the behavior 
of a typical solution. A typical solution is the superposition of all four solutions so 
that the largest exponent dominates. 



s 




Fig. 4: The dark line shows the behavior of a "small solution" as we change 
the argument of z. The solution is controlled by the smallest exponent in the 
darkened angular sector. The solution is still controlled by the same exponent as 
we move away from the darkened region, until the exponent grows to a maximum. 
Subsequently, it behaves as the generic solution, see fig. 3. 

There is a simple, powerful fact: if we pick a ray towards |z| — )• oo away from the 
boundaries of the Wj, we can identify a unique (up to rescaling) exact solution by the 
boundary condition which requires the fastest possible decay at infinity. This small solution 
Si is the same for all rays inside Wj. More generally, one can identify uniquely A;-dimensional 
subspaces of solutions whose growth is controlled by the k smallest exponentials. That 
way, it is easy to see that Si will still behave as exp^i if continued beyond Wj, through 
Wi-\ and Wi-i or through VFi+i and as 5"^ grows all the way to be the largest 

exponent see fig. 4 . Beyond that, the behavior of Sj is unconstrained. 

One can generically take the set (sj, Sj+i, Si+2, Si+a) to be a neat basis of linearly inde- 
pendent solutions. This implies that can be expressed in terms of a linear combination 
of the previous four solutions. In the sector where s^j^i is the smallest solution, Sj and Sj-1-4 
are the two largest solutions. Thus we can compare them. We choose to normalize so 
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that it is proportional to —Si in this sector, up to corrections by smaller solutions. Thus 
we can write 

Si+4: + Si = ai_(_iSi_(_i + bi+2Si+2 + Ci+3Si_|_3. 

With this choice of normalization we have that det(st, s^+i, Si+a) = 1- We should 
impose the periodicity constraint = fXiSi. We can relate the /i^ with the formal 

monodromies, i.e. the monodromies of the approximate asymptotic solutions. The formal 
monodromy receives contributions from both m, m and the coefficient in front of 1/z in 
the expansion of P(z)^/'*, if present. The simplest example is the pentagon. There are no 
formal monodromies, and Si + S2 + S3 + S4 + S5 = is the unique way to solve all the 
constraints. 

The example of interest for us is the hexagon. In this case our relation requires si, S5, 
■§95 '§13 = /^I'Si to grow with the same exponent Si{z), but the comparison is made after 
going around the z plane twice, hence nf — exp (5*1 (e^^*^) — Si{z)), and fxi agrees up to 
a sign with the formal monodromy. The formal monodromies come from m, rfi ( |2.21j ) only. 
We can set without loss of generality P{z) = z'^ — U, and then P(z)^/'* = 2;^/^ — ^^72 + ■ ■ • 
has an expansion in half integer powers of z~^ and there is no logarithm!! in w. Hence 
we get formal monodromies S7 = fxsi, sg = /U~^S2, sg = fxss, sio = /U~^S4, etc. Here 
(U = ± exp 27rz(m — m). In (3, 1) signature, /U is a phase, while in (2, 2) signature ^ is a real 
number. 

The coefficients hi will play a particularly important role in what follows, since 
for n = 6 they contain the non-trivial Stokes information. Consider the determinant 
det(si, Si+i, Si+3, Si+4) = —bi+2- For the hexagon that is the same as 

bi+2 = -/U/U""^ det(Si+6, Si+3, Si+4) = bi+5 

It turns out that everything can be expressed in terms of fi, bi, 62, ^3, subject to the 
relationffl 

6162&3 = bi + b2 + b3 + ^ + h+3 = bi (2.27) 



^ More explicitly, in cases where A'^ is a multiple of four, we will get that generically we 
have dw = {z^^^ + • • • + rh/z + • • ■)dz which implies that w contains a term going like w ~ 
• • • + m log z + ■ ■ ■. Such a term also contributes to the formal monodromy. 

^ We can change bi — >■ —bi and /x — > by changing the signs of some of the Si. Thus these 
two possibilities describe the same physical situation. In this paper we concentrate in the regime 
where the sign of all biS are the same, which we take to be positive. Then the two signs of /i 
describe two different physical situations. 

^ The /U = 1 limit of this equation was obtained in the study of the anharmonic oscillator 
in [|14|. 
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The case of the hexagon in AdSs is particularly important. Then, the solutions at 
even and odd sectors do not mix, and satisfy si — S3 + S5 =0 and S2 — S4 + sq = 0, so that 
hi = 1 identically, and /U = — 1, as it should. For the AdS4^ case we have that /U = ±1. 

To conclude this subsection, we would like to relate the bi to the space-time cross- 
ratios for the case of the hexagon. If we are given a solution qi of the vector problem, 
there is a very simple way to extract the position of the cusps. The bilinear sJV^ Si^^i 
gives the vector solution which decays the fastest in the sector Vi, hence the contraction 
vt = r^Si+i gives the position of the cusp. The inner product is 

Vi ■ Vj = (Qi ■ Qj)s'[T^Si+isjT^Sj+i = sfV^ Si+isJViSj+i = det(si, Si+i, Sj, Sj+i). 

Hence the cross-ratios are 

™2 ,yi2 1 ^2 ^2 1 /yi2 ™2 1 

In conclusion, we have shown that the Stokes data at infinity is given by 61, 62, ^3 
and fx with the constraint ( 2.27| ). The spacetimes cross ratios are related to the bi via 

(Pi). 

3. From Hitchin to TBA 

If we introduce now the spectral parameter, C, in the connection, we can investigate 

the behavior of the deformed solution Y[Q. As the conserved polynomial becomes simply 

C~^-P(^), the deformed solution will have the same number of cusps at infinity, associated 

with directions yj[C] and cross-ratios Uj[Q (or bj[Q) which are meromorphic in C away 

from C = 0, 00. The only slight subtlety is that the location of the sectors Wj[(] and Vj[C,] 

varies as the phase of C varies. The sectors Wj [Q are defined by 

r>n 27rj 4 ^ Stt 2nj 4 ^ 

C '■ 1 arg Q < arg z < \ arg Q . 

n n n n n n 

Hence the small sections Sj [Q , the coefficients bj [Q and the cross-ratios Uj [(] do not quite 
come back to themselves as ^ — ?■ e^^*C, but rather undergo a shift in the index j ^ j + 4, 
bj[Ce'^'^'^] = bj+4[C]. Due to the Z4 automorphism of A, $, there is actually a symmetry 
bj [iQ = [Q ■ We will not make use of this symmetry for now. Note that, for the hexagon 



case, the relation (|2.27| ) that we had for C = 1 remains valid for bi [Q , with 11 independent of 
(. Throughout this section we will hold U and ^ fixed and we will only vary (. Of course, 
we are only interested in the values at ^ = 1. We will determine these values by writing an 
equation for the cross ratios as a function of ( and then we will set C = 1- This equation 
comes from demanding that the functions are analytic plus the statement that they have a 
known approximate expression (displaying Stokes phenomenon) near C ~ 0, 00, which can 
be computed via the WKB approximation. 
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Fig. 5: This figure siiows tiie sectors in the z plane for the case of n = 6. They 
could be the sectors Wj or Vj. In (a) we see the various sectors for a given value 
of the phase of C- (b) as we change the phase of C by C — >■ the sectors rotate 
by an angle ^ipjZ. In (c) we see the result of a full rotation C, — >■ ^e^^"^. We come 
back to (a) up to a relabeling of sectors. 

We aim to set up a Riemann-Hilbert problem, and derive a useful set of integral 
equations from it. The equations take the appearance of TBA equations. Although the 
resemblance is purely formal, it is rather striking. The TBA equations do not just compute 
the cross-ratios. We will show how the regularized area of the minimal surface, which can 
be written as a suitably regularized integral of Tr$$, takes the form of the free energy for 
the TBA equations. 

For the sake of simplicity, we will give a direct derivation of the TBA equations for 
the case at hand. The reader should be aware that this is a special case of a much broader 
result. The mondromy data (such as the cross-ratios) of the auxiliary connection for any 
Hitchin system is always computed by a TBA-like set of equations. Physically, this fact 
follows from a relation between Hitchin systems and M = 2 gauge theories pB| , pU| . In 
the special case of SU (2) Hitchin systems, the result has been demonstrated directly, and 
a crucial ingredient for the TBA equations (the spectrum) has been determined through 
a simple WKB analysis ||2^. For more general higher rank systems, the simple WKB 
analysis is expected to fail, and be replaced by something more complicated, and currently 
unknown. Luckily for us, the special Hitchin systems with a simple spectral curve = 
P{z), such as the one relevant for the minimal surfaces in AdS^^, are still amenable to a 
simple WKB analysis, although an extra layer of analysis is required. The relation between 
the regularized integral of Tr$$ and the TBA free energy appears to be valid for general 
Hitchin systems as well, and was not noticed in previous works. 
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3.1. WKB approximation o/6j[C] 

Near ^ = 0, oo the behavior of the connection is clearly dominated by $ and $ 
respectively, and we can use a WKB approximation to evaluate the behavior of the cross- 
ratios. The basic idea behind the WKB approximation is to go to a (complex) gauge where, 
say, both $ and A have been diagonalized. If one ignores the higher order corrections as 
C — > 0, the z-th component of a flat section should behave like exp ^ J Xadz., where Xa are 
the eigenvalues of The higher order corrections mix the different components of the flat 
section, and disrupt this simple exponential behavior. The component with the highest 
exponential growth (largest real part of ^) will in general contaminate all the others, 
but its leading exponential growth will be unaffected. The simple WKB approximation 
computes the small C, exponential behavior of a flat section transported along a path for 
which the real part of ^ for a given eigenvalue remains the largest. This is called a "WKB 
path" : the condition guarantees that the flat section grows with the largest exponent and 
is well described by the WKB approximation. 

In the case at hand, we want to compare the "small" flat vector-valued sections 
sfr^Si_|_i at the rays ri between the sectors Wi and VFi+i. At sufficiently large 2;, the 
small section is, by definition, the flat section which grows the fastest along a line which 
goes towards smaller z along that ray r^. If we can flnd a WKB path which flows all 
the way to some other ray Vj at large 2;, we can evaluate the small C, asymptotics of the 
inner product det(si, Sj+i, Sj, s-,+1) as exp ^ Jj;' xdz. pij are some reference points near 
2; = 00 in Tij. Varying pi, pj changes the (arbitrary) choice of normalization for the small 
sections. Remember that the eigenvalue x for the vector- valued connection will be the 
sum of two consecutive roots P(2;)^/^. If we have appropriate WKB paths available for all 
the pieces of a cross-ratio, and we can evaluate the asymptotic behavior of each piece and 
combine them. The spurious dependence on the pi drops off from the cross-ratio, and the 
various contour integrals join into a single integral of P{zy/^ along an appropriate contour 
on the Riemann surface x"^ = P{z). 

The eigenvalues Xa for the vector valued $ are of the general form 
0, 0. The real part of C,~^xdz is the largest if the phase of C,~^xdz sits inside an open 
sector of width 7r/2 centered around the positive real axis. Beyond that sector, either 
C,~^ {—ix)dz or C,~^{ix)dz will have a larger real part. This is the condition which must 
be satisfled by a WKB path. It is particularly useful to consider WKB paths of steepest 
ascent, such that the phase (p of xdz is constant. Such a path satisfles the WKB condition 
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for all values of ( inside a sector of width 7r/2 centered around the ray of phase cp in the ( 
plane. Correspondingly, any estimate of a cross-ratio asymptotics done with WKB paths 
of constant phase (p will hold along all rays inside the corresponding sector in the ( plane. 
Notice that the WKB condition is satisfied in the sector both with respect to C,~^xdz and 
to C,xdz, so the same calculation controls the asymptotics both at small and at large |C| 
for a given phase of C- 

The WKB paths of constant phase </? are straight lines in the w plane. At sufficiently 
large z, away from the zeros of P{z), there are always WKB paths joining rays and ri_(_2 5 
which give estimate of the reference inner products {si A A Si_|_2 A Si+3). The straight 
lines in the w plane can be grouped in continuous families which asymptote to the same 
rays at infinity. Different families are separated by the special lines which end on the zeros 
of P{z). It is useful to draw a picture in the z plane, with one WKB path from each family 
for a given choice of ip (and all possible choices of a fourth root of P{z)). 




Fig. 6: This figures shows WKB lines in the z plane with a constant phase. The 
doted lines end at the zeros and they divide various regions. WKB lines that cross 
at 90 degrees correspond to two different eigenvalues that differ by an i factor. 
SoUd fines show fines that we wifi use to evaluate 63 and 61 in this case. The figure 
on the right shows the effect of increasing the value of if. Note that we have WKB 
lines joining cusps 12 with 45 and 23 with 56 but we do not have WKB line joining 
16 with 34. 

We see that only two out of three opposite pairs of rays are joined by a WKB path, 
which allows the estimate of {si A Sj+i A Si+3 A 5^+4). See fig. 6. As one varies the value 
of continuously, the configuration of paths of constant phase also varies continuously. 
Jumps happen at special values of where the special WKB paths ending on a zero of 



19 



Fig. 7: These are the same WKB lines as in the previous figure but for increasing 
values of ^p. The figure on the right represents a critical value of the phase where 
the line that was joining the cusp 23 with 61 no longer exits. Instead we have 
WKB lines that end at zeros and a new WKB line that joins the two zeros. In 
this case we can only use the WKB approximation to compute the cross ratio h\ 
involving cusps 12 with 45. 



Fig. 8: This figure shows the WKB lines for a phase which is bigger than the 
critical phase. In this case we have WKB lines joining the cusps 12 with 45 and 
also joining 16 with 34. In this case we can compute 62 and 63 via the WKB 
approximation. 

P{z) collide with the other zero, see fig. 7. Clearly, the special values of (p coincide with 
the phase of the integral of xdz from a zero to the other along that path. Since x is the 
sum of two eigenvalues of the spinor connection, we have x = ^'^^ [z^ — UY^^, see ( |2.20|) . 
We can rescale z — ?■ U^^'^t to get xdz = U^/'^{1 — t^Y^'^dt. The path runs on real t between 
t = —1 and t = 1, the special values (pi are the phases of the four roots of t/"^/^. At each (pi 
a single pair of opposite rays (which we can identify without loss of generality with and 
r^-i-s) are joined by a WKB path of constant phase. The path exists in the whole range 
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(Pi-\-i < if < ipi-1, i.e. {pi — TT/2 < (fi < Lpi + iT/2. Hence the product {si A s^+i A s^+a A 5^+4) 
has a simple WKB estimate for a range of 

(pi - 3n/4 < argC < (fii + Sn/A. 

We can use this information to estimate the behavior of bi[Q as C — > and C — > 00 
inside this sector. We find exp Zi/( aX smaU |^| and exp at large |C| on the whole sector 
(Pi~3TT/4 < argC < Lpi + Srr/A for an appropriate function Zi[U]. Patiently assembling the 
WKB phase integrals for all the pieces of the cross-ratio, Zi is 



3r(|, ^ '^-^ 

The different Zj corresponds to the four roots of t/^/^. The specific root of C/^/^ is exactly 
the one with phase ipj , hence bj is exponentially large in the half plane (pj — n /2 < ip < 
(fj + tt/2. The asymptotic behavior of bj changes beyond the lines at (fj ± 37r/4. We can 
check this explicitly using the relation ( p.27|) . For example, we can write 

bj = + ^"^ + bj-i + bj-2)/{bj-ibj-2 - 1). 
As we cross the ray at <f>j + 3tc/4, bj-i and bj-2 switch dominance, and the asymptotic 
value of the ratio passes from bj ~ 6~^2 ~ ^^P[^j/C] to bj ~ ^J-i ~ 6xp[— ~ 
exp[—iZj /(]. On the other side, beyond the ray at (pj — 2>tx/A we get bj ~ b~^^ ~ 
exp[-Zj+i/C] ~ exp[iZj-/C]- 



bi 



Fig. 9: We are displaying the C, plane for a particular value of the phase of 
U, chosen so that along the positive real axis, b\[Q becomes largest. Along the 
positive imaginary axis, h^[C\ is the largest, and so on. Due to to the fact that 
the bj[C,e^'"'^] = 6j+4[C] = ^j+i[C] we have a branch cut, denoted by the doted line 
where we switch the indexing of the bj. We have indicated what each bj changes 
into. As we change the phase of U the lines where bi are the biggest rotate rigidly. 
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This is all the information we will need from the WKB analysis. 



3.2. A Riemann-Hilbert problem 

In this subsection we set up a Riemann-Hilbert problem. For this we need to combine 
the fact that the bj[Q are meromorphic in ( (except for C = 0? oo) with the fact that they 
have different asymptotic behaviors as ^ ~^ 0, oo depending on the phase of C- In other 
words, the bi[Q display Stokes phenomenon, but now in (. 

In order to set up this Riemann-Hilbert problem, we would like to identify some 
functions Xa of the bi whose asymptotic behavior at small and large ( is uniform on the 
whole C plane. More concretely, we would like to define the Xa in such a way that they 
behave as Xa ~ e'^"/^ as C — )■ 0, for all phases of C- Here Za denotes the four eigenvalues 
of the connection which go like Z, —iZ, —Z, iZ. Of course, due to the Stokes phenomenon 
in bi, the price to pay is that there are discontinuities in the Xa[C] along some rays in the 
C plane. The Riemann-Hilbert problem is captured by the data of such discontinuities. As 
we will show explicitly below, such discontinuities can be extracted from the various WKB 
approximations discussed above. One may set up many slightly different Riemann-Hilbert 
problems, based on slightly different definitions of the Xa. There is a canonical choice 
though, inspired by the connection to A/" = 2 theories in four dimensions , which gives 
discontinuities well suited to formulate the final integral equations we are aiming for. 

It is instructive to first consider a non-ideal set of functions T^, allowing for disconti- 
nuities at the rays ipi in order to keep the uniform asymptotics ~ expZa/C uniformly 
true everywhere. For instance, let us define 

Ti[C]=6i[C] for -^<argC-(/Pi<^ 

Q (3.2) 
Ti[C] = l/62[C] for ^<argC-^i<y 

where ipi is the phase of Zi. We have similar expressions for T2 where we change bj — ?■ bj^i 
and fi^f2 = fi— 7r/2 in the above expression. We then define Ta+2 = '^a^- Now, we 
would like to show that has uniform asymptotics. Let us start with Ti. For simplicity, 
let us set the phase of Zi to zero, we will restore it later. In the upper half ( plane 63 is very 
large and we can use ( p. 271) to write 6162 — 1 = ^3 ^(^i+^2 + A* + 1/a*)- As 63 is exponentially 
larger at small ( than the factor in parenthesis around the boundary arg^ = ipi + n/2 
(see fig. 9), we see that 62 ~ l/^i- This shows that I/62 ~ exp[Z'i/(^] near this line. But 
since the asymptotics of 62 is already given by the WKB approximation, we see that this 
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continues to be true throughout the second sector in (|3.2|). The same is true for all the 
T 

-■- a- 

We aim to derive integral equations for the logarithms log T„ from their discon- 
tinuities. At this stage, the discontinuities are rather ugly-looking. For example, at 
argC = -f 7r/2, Ti goes from bi to 63 ^, and the ratio of the two values is a complicated 
function. This can be easily computed by using the relation ( |2.27| ) to write 

Ti(c+)/Ti(r ) = 61/62"' = 6162 - 1 + 63-1(61 + b2 + ^i + ^i-'). 

Here denotes arguments for C which are slightly bigger or smaller than the discontinuity 
line. It cannot quite be written as a combination of which is continuous on this line. It 
is always convenient to write the discontinuities in terms of functions where are continuous 
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Fig. 10: Here we see the explicit definitions of Xi and in each angular sector 
of the C, plane. (We have set the phase of Z\ to zero for simplicity). Of course the 
bi are discontinuous only at the branch cut. The fact that we define Xa differently 
in different sectors is the source of the discontinuities. 



We now define new variables Xa which have simpler discontinuities, although more 
discontinuity lines. Let us define Xa{C) to coincide with ba only for — 7r/4 < arg(^ < 
fa + T^/^:- To define it in the remaining regions we define XaXa-2 = 1 and also require 
XaXa-i = baba-1 — 1 for < arg ^ < (^a-i- The important property of XaXa-i is that its 
analytic continuation has a simple WKB expression in a sector of angle slightly bigger than 
TT. This defines each Xa on all the quadrants. In fig. 10 we see the resulting definition of 
Xi, Xs. Compared to the variables Ta we have now introduced extra discontinuities. The 
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Fig. 11: We see how Xi and X3 change as we cross the lines of discontinuity. 
We have defined = (1 + fi/Xa){l + Aa,a-i = (f + x^x^_J - These are 

defined so that, for example, across the top vertical line we have = A^^ 
where denotes the value of just after crossing the line (and just before) 
if we move in an anti-clockwise fashion. 



correction is mild, but important. There are now eight lines of discontinuity. At ipa, '^a is 
continuous. For Xa-\ we have the definitions X^-x = 6a-i — ^/ba and A'a+i = ^>a+i — 
on the two sides of the (fa line. We find 

where we used (|2.27| ). Hence the discontinuity of \ogXa-i is a function of the continuous 
variable Xa- log(l + ijlX~^){1 + ijl~^X~'^). At + 7r/4, XaXa-i is continuous by design. 
For Xa, we need to compare the definitions Xa{C~) = ha and Xa{C,'^) = ha — 1/ha-i- We 
find 



haha-1 

Hence the discontinuity of log X IS 3j function of the continuous variable XaXa-i only: 
— log(l -f -^-^ — -). The resulting discontinuities are summarized in fig. 11. 

We have our final Riemann-Hilbert problem: two functions logA'afC], with uniform 
asymptotic behavior at large and small and discontinuities which are functions of the 
Xa themselves. The discontinuities take the canonical form expected from the connection 
to A/" = 2 theories in four dimensions. Indeed, both the asymptotic behavior and the 
discontinuities of the Xa{C) are identical to the ones for the cross-ratios of an SU(2) 
Hitchin system, associated with the = 4 Argyres-Douglas theory in the notation of 
|2y[ |. The spectral curve for that system takes the form = + U, which is identical 
to the curve for the hexagon if one takes z ^ x. The transformation also relates the 
canonical differentials xdz, as xdz + zdx is a total derivative. This relation guarantees 
the coincidence of the periods Zi, but it does not imply the identity of the discontinuities. 
That is a non-trivial, and rather unexpected fact. 
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3.3. An integral equation 

Let us first assume that the functions Xa[Q have no poles or zeros in their region 
of definition. Along various lines the will have discontinuities that follow in a simple 
way from the above definitions. In fig. 11 we have summarized the information about the 
discontinuities both for Xi and X^. 

Functions with such discontinuities can be obtained as the solution of integral equa- 
tions of the form 

+ ^n^,K„ .log ((1 + + _!_)) + + 

(3.3) 

The Kernels are making sure that when ( crosses the contour we pick a contribution which 
is precisely the one we expect for the corresponding discontinuity. The integers Uab, triab 
can be read off from figure fig. 11. Namely, for a = 1 or a = 3 they are the exponents of 
the Ys in fig. 11. They have the values 

Uaa = na,a-l = "1 na,a-2 = na,a-3 = 1 

(3.4) 

rUaa = -1 TOa,a-l = "1 W„,a-2 = 1 "^a,a-3 = 1 

This equation can be further simplified by using the fact that A'^-ifiC] = ^a[C] which 
is due to the Z4 symmetry of our Hitchin problem. This implies that the values of Xa that 
are appearing along the rays are all the same function. Let us set the phase of Zi to 
zero for simplicity. We can then introduce two functions e and e which parameterize the 
values of the Xa and XaXa-i at all the rays. Namely 

e<^^ =Xi[C = e^] , e^'(^) = XiXs[C = e^'/^e^] (3.5) 

The functions Xi and X^ on the full ( plane can be determined by the two functions of a 
real variable e and e through the right hand side of ( |3.3| ) . The integral equations for e and 
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(3.6) 



e follow then simply by evaluating the left hand side along the rays used to define e and e. 
We thus get 

+ ^ / Ufl .., log(l + ^e-)(l + — ) 

ZTT y cosh^y — 0') fx 

m =2V2\Z\cosh9 + - ! de' J- log(l + e-^") + 

TV J cosh^y — 6') 

These are completely explicit integral equations for two functions. The integrals run 
between minus infinity and infinity. The functions e, e inside the integrals are evaluated 
at 9' . These turn out to coincide with the TBA equations of the As (or Z4 symmetric) 
integrable theory introduced in and further studied in |]^. This theory has two 
particles of mass m and one bound state with mass very similar equation was 

obtained in the study of an x"^ anharmonic oscillator in . Their equation can be obtained 
by the replacemnt \Z\ cosh^ — )■ which is related to the UV, or conformal, limit of the 
As theory. 

Once we find a solution to ( |3.6| ) we can then write the cross ratios, bi, which we 
can define through the right hand side of (|3.3|) and from these we can construct the Ui. 
Denoting Z = \Z\e''^ then note that ha[Z, C = 1] = ba[\Z\X = e-'^P]. We can pick a definite 
octant, say — 7r/4 < ip <Q. We can then compute all three bi on this octant from the two 
expressions 

log(— -l)=logA'iA'3 = 2V2cos(^|Z| + - / dO' I . log(l + e-^") + 
U2 Tt J cosh(6' + lip) 



+ — / d6' — — ^— — ^log(l + e + 



TV 



cosh 2(6*' + i(p) n 



log 61 =log;^i =2|Z|cosv.+ ^ / de' ^"fi;^ + ^I log(l + e-') + 



(3.7) 



TT J cosh.2{iip + 6' 

2tt J cosh(z(/? + U'j /i 
Where ip = f + n /4. From these and (|2.27|) we can determine 

1 



b2 



biU2 

bi+bs + fi + j^ 
bibs - 1 
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(3.^ 



Once we have defined the bi on this octant, we can use the symmetries of the problem to 
map them to other octants. In fact this single octant is a complete solution, since we can 
always permute the ui so that we are on this octant. 

Notice that the equation we wrote is almost indifferent to the signature. The only 
difference is that in (3,1) and (1,3) signature /i is a phase, while in (2,2) signature fx 
is real. The main question one may ask is if the condition on the absence of zeros and 
poles in the Xa is self-consistent. There is a regime where the self-consistency is manifest, 
namely ^ 1. In that case the WKB approximation is valid essentially for all and 



Xa ~ exp 



C 



+ ZaC 



are such that the logarithms on the right hand side of ( |3.3| ) are 
exponentially small along the contours of integration. In the set up analyzed in this 
corresponded to a certain "large radius" limit. Thus the problem becomes simple for 
1^1 ^ 1 and it is a good starting point for a numerical solution of the equations. It 
corresponds to the low temperature limit of the TBA system. 

Let us see what large \U\, or large \Z\, corresponds to in our case. Let us assume that 
we are in a regime where bi and 63 are both large and have a good WKB approximation. 
Writing Z = \Z\e^'^, this corresponds to —tv/2 < ip < 0. We then have the expressions 

TT 

log 61 ~ 21^1 cos(p , log 63 = 21^1 cos{ip + — ) 

11 ^ ui bs 

02~^ + 7 ^wi+ws^l, — = 7- (3.9) 

61 63 Us bi ^ ^ 

1 

U2 ~ -r-r- < 1 

bibs 

We see that if the phase (p is generic then we have (mi,W2,'U3) = (1,0,0) or (0,0, 1). 
On the other hand if the phase (f> approaches — 7r/4 as \U\ — > 00, then we can have a finite 
value of ui/us and we get a segment (1 — W3, 0, W3), < W3 < 1 which joints the two points 
mentioned above. In other angular sectors we get other segments with Ui — ?■ w^+i. 

Thus, the large \U\, or large \Z\, limit corresponds to cross ratios that are on this 
triangle. We can move away from this triangle by lowering the value of \U\. U, together 
with /i span the three dimensional space of cross ratios Ui. Within this space there will be 
solutions that can be interpreted as living in various signatures. In the next subsection we 
describe more precisely the region in the space of cross ratios which is spanned by polygons 
living in space with various signatures. 
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3.4- Kinematics, signature and the space of cross ratios 

In this section we discuss the kinematics relevant for the hexagonal Wilson loop (or 
scattering of six particles), in several signatures. We consider an hexagon with cusps at 
points Xi in four dimensions, such that the distance between two consecutive points is 
light-like. One can construct three independent cross-ratios 

^1 - :3-;2-' «2 - -2—2-' ^3 - -2—2- (^-1") 
•*'14**'36 •^25-^14 •^36-^25 

we would like to understand the range of such cross-ratios in different signatures, namely, 
(3, 1), (2, 2) and (1, 3). It is convenient go to a conformal frame where we have sent three 
points, for instance a;6,i,5, to infinity. Furthermore, we set the position of Xs as the origin. 
To be more specific, we denote four dimensional coordinates by a; = {x~^ ,x~ ,x±) and 
choose 

X3 =(0,0,0), X2 = (-1,^2 X4^{x^,-l,q) 

(3.11) 

a^i =3^2 + (0,A,0), a;5 = a;4 + (A,0,0), = (A, A, 0) + .... 

Since X2 and X4 are light- like vectors, we require X2 = —p^, x'^ = —q^. Here, p and q 
are two vectors living in the transverse space. This transverse space has various signatures. 
They can be used to write the cross ratios. 

1 1 + g - 2g.p 1 , . 

«i = --. ^, «2 = ^— W3 = Z2. (3-12) 

1-q (1 - g^)(l -p ) 1-p 

The regular hexagon embedded in AdSs corresponds to q — p — 0, with Ui = 1. We are 
interested in solutions continuously connected with the regular hexagon, for which Ui > 0. 

A negative sign for some Ui means that some of the distances between cusps are becoming 
timelike. In this paper we will discuss only the case where the Uj > where distances are 
not yet timelike. 

Let us discuss in more detail the physical meaning of the large \U\ lines, for example 
U2 — 0, u\ + — 1. The second condition sets q p — 1. Since, Ui are positive, we see 
that Ml, Us are less than one. This means that q and p are negative. The condition 
U2 = requires p ■ q =1, which implies that p and q are collinear and point in opposite 
directions. We see then that the four dimensional vectors X2 and x^ are collinear as well. 
We conclude that the large \U\ regime is actually the collinear limit. In fact, the fraction 
^ — ^fcn^ = where the ki^i+i = Xi+i — Xi are the spacetime momenta. Thus, as we 
move along the segment (1 —1/3, 0,143) we change this fraction which is a free parameter in 
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Fig. 12: This is the three dimensional space of cross ratios parameterized by 
(til, U2U3). The lines belonging to the edges of the triangle correspond to the large 
\U\ limit. These edges correspond to various collinear limits. The vertices of the 
triangles correspond to soft limits. 



the collinear limit. Note that for z to be different from zero or infinity we need to fine-tune 
the phase of as we move to large values of \U\. However, in this regime the area, A free 
is not sensitive to this phase. In fact, we will see that our remainder function remains 
constant in this limit. Taking the large \U\ limit along a generic direction gives us a simple 
soft limit which simply eliminates one of the momenta. Let's say that u^^ui — )■ in this 
limit, then we conclude that q — > 00, which after a rescaling , is making X2 ~ xs- Thus, 
we are sending the momentum /C23 — 0. 

Note that the expressions ( |3.12|) are valid for any signature under consideration. Of 
course, the sign of depends on the signature. For instance, for (3, 1) signature, > 
and hence Ui > 1, while for (1,3) signature, < 1 and Ui < 1. However, the converse is 
not true. The reason is simple, with (2, 2) signature we can have either sign for q. 

To get more information we look at the sign of p^q'^ — [p-q)'^, which is positive in 
(3,1) and (1,3), but is negative in (2,2). This can be translated into a condition on fx. 
Combining (|2.27|) with (|3.10|) we can write n as 



Using the above formulas we can write 



1 ~ Ui - U2 - Us 



-1\2 



y/UiU2U3 



U1U3 /^2-2 
[P Q 



(3.13) 



(3.14) 
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Fig. 13: Here we see the AdS4 surface corresponding to ^ = ±1 in the space of 
cross ratios parametrized by {ui, U2,us). Points on this surface corresponds to cross 
ratios that arise when we have a polygonal contour in a three dimensional space. 
The surface pinches at iti = 1. The region inside the "bag" in the Ui < 1 region 
correspond to configurations in (1,3) signature. The region outside corresponds 
to surfaces in (2, 2) signature. Finally the region inside the surface for Ui > 1 
corresponds to ordinary (3, 1) signature. We show the same surface from two 
different viewing directions. 

We conclude that is a phase for (3,1) and (1,3) but fx is real for (2,2). Thus the 
boundary of the (2, 2) signature lies at = ±1. These values coincide with the ones 
corresponding to surfaces that can be embedded in AdS^. This is not surprising, we need 
to make one spacetime coordinate vanish if we are going to change its signature. We plot 
this surface in fig. 13. The surface pinches aX ui = U2 = us = 1. This corresponds to the 
regular hexagon embedded in AdSs, all three signatures meet here. The AdS^ surface has 
two sheets emanating from the large \U\ triangle. One sheet has /i = — 1, it starts at the 
large \U\ triangle, goes through the pinch at = 1 and then continues through to Ui > 1. 
The other has fi = 1 and it stays at < 1. It has the topology of a disk and it contains 
an interesting polygon with ui = U2 = us = \- This is another kind of regular hexagon 
which has a simple picture in (2, 2) signature and it is discussed further in appendix B. 
These two branches together define the full AdS^ surface. The region inside the surface 
for Ui < 1 has (1, 3) signature, the region outside has (2, 2) signature and the region inside 
for Ui > 1 corresponds to ordinary (3, 1) signature. 

It is interesting to look at the surfaces of constant which end on the large radius 
triangle. We have already discussed the fx = ±1 surfaces. If is a phase the surface has 
the topology of a disk which sits in the (1, 3) signature region (one space three times). We 
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expect that as we solve the integral equation for decreasing values of starting from the 
large radius region, we will get a nice map from the U plane covering the disk once. It 
turns out that U = corresponds to the line ui = U2 = u^. We move along the line by 
changing fi. In section five we find the exact solution corresponding to this case. 

If we take fi to be real and negative, the surface has the topology of a cylinder in the 
region Ui > 0, and has some pieces in the regions with some Ui negative. We expect that 
as one moves towards small \U\, it will cover the cylinder, and then at some point one 
ba will become negative, making two Ui negative. In order to reach the (3, 1) signature 
region, a possible strategy is to solve the integral equation starting from the large radius 
region, with fj, real and slightly smaller than —1, follow it until we get to the Ui > 1 region, 
and then cross into the 3, 1 signature by changing n back to a phase. Some interesting 
phenomena can occur in the process, if the functions Xa[Q acquire zeroes or poles. This is 
a common (even desirable at times) occurrence in solving TBA-like equations p2| , |33[1 . In 
appendix C we discuss this a bit further. 



4. Area and TBA free energy 

In this section we will show how to compute a certain regularized area as the free 
energy of the TBA system for the hexagon case. 

The area is given by 2 J e" which is divergent. This should be regularized in a physical 
way. For the time being let us just define a regularized version by subtracting the behavior 



of e" far away. We have seen that a (|2.19| ) goes to zero at infinity. This means that for large 
e" ~ Furthermore, in the case of the hexagon, we also have that ~ \z\ at 

infinity. Thus we find it convenient to introduce the following definitions 

-^reg -^periods ~^ -^free 2 d (6 |^|) 
Aperrods =2 j (f z{\P {z)\^ - \z\) = \Z\^ (4.1) 
Afree =1 j ^^^(e" - l^l^/^) 

where we have evaluated explicitly Ap^riods in terms of ( |3.1j ) . Note that this is a convergent 
integral in our case since P{z) = z"^ — U . 

It is convenient to write A^gg in terms of the "Higgs" fields of the Hitchin system 
using 

Tr$^$^ = e° - v^vie-"^ = 2e" - dda 
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Hence, up to a total derivative which integrates to nN , we can define (the trace is in the 
spinor representation) 



A 



free 



j d^z[Tr^,^, - 2\z\]d'^z 



The formal expression J Tr^^^z has a nice interpretation in the context of the Hitchin 
system. Recall there is a natural Poisson bracket on functionals of A, $ 



SU dV su sv _ 

s¥:,6¥;^^6a:,6a;.~^'^'' 



(4.2) 



and the fiow generated by / Tr^^^z rotates $2 and $2 in opposite directions as 5$2 = 



1$ = 
2^^ 



The Poisson bracket is gauge invariant, and it actually descends to a nice Poisson 
bracket on functionals of the fiat sections of the spectral connection, such as the bi. On 



such functionals, it reduces to 



V} = i J (fzTr 



6U 6V 



- [z ^ z) 



(4.3) 



It is straightforward, but tedious, to compute the Poisson bracket {bi, bi-i} = bibi-i — 
1 along the lines of appendix B of [^. There is a self-consistency check, one can see that 
{bi [Q , [C] } had to be some polynomial in the bi [Q with constant coefficients independent 
of C. The right hand side bibi-i — 1 is the only possibility consistent with the asymptotics 
of bi at small ( along all the rays in the ( plane. Also, {6^, n} = 0. Another consistency 
check is that this Poisson bracket is consistent with the relation (|2.27| ) . Taking the Poisson 
bracket of the two sides of ( |2.27| ) with bi we get 

{61, 616263} = 6162(6163 - 1) - 61(6162 - 1)63 = 6163 - 6162 = {61, 61 + 62 + 63 + ;U + IJ~'^}. 

Finally, the neatest result: {Xa, Xa-i} = XaXa-i identically in all sectors. In particu- 
lar, the discontinuities in the Xa[C] preserve the Poisson bracket, and the corresponding 
symplectic form 

u = dlog A'ldlog ^"3 = dlog Afi(ilog(A'iA'3) . 

The asymptotic behavior of A'a at ^ = 0, 00 shows that u holomorphic in C everywhere, 
and hence is ( independent. 
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A rotation of and in opposite directions can be absorbed by a phase rotation 
of C- Can we identify a function F which can generate this phase rotation directly? It 
should satisfy the equation 

dF='-[d log A'l Cdc log Xs-d log C^c log X^] . (4.4) 

F , as will be ( independent. This allows us to compute it at, say, very small From 
the integral equation, we can extract the behavior of log Xa at small ( 

logXa^^+Ka + C {Za + K) (4.5) 

Here Ka and K'^ are simply the value at C = and the first derivative at C = of the 
convolution terms in the integral equation. Then 

and dF can be written as dF = dFi + dF2 where 

dFi = ^d{ZiZ3 - Z1Z3) = d\Zf (4.6) 

where we used that Z\ = Z and Z3 = iZ. This is the dominant contribution in the large 
\Z\ limit. The second term is 

dF2='-d{Z^K'^-Z^K'^) (4.7) 
F-2 can be rearranged a little, to the form 

(4.8) 

A similar calculation at large C produces a similar result, but with ^ replaced by ZaC- 
The two expressions must be identical, and we can take the average 

Notice that F2 goes to zero in the large \Z\ limit. In this limit the dominant piece is 
Fi , which reproduces ^periods • Note that there could have been a fj, dependent integration 
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constant in F. However, this is not the case because we know from the properties of the 
solutions that A free as define in (|4.1|) goes to a constant for large \Z\. Thus the conclusion 
is that, up to a constant, we have that Af^ee = -^2- Finally, let us write Af^ee in terms of 
the solutions of the integral equation ( p77|) 



^free = — / rf^2|Z| cosh ^ log (1 + e"'/x) (1 + — ) + 2v^|Z| cosh ^ log (1 + c""') (4.10) 

This has precisely the form of the free energy (up to an overall sign) of the TBA system 
in (pD. 

In order to get the full physical answer of the problem we need to introduce a more 
physical regulator. The physical regulator consists in cutting off the integral when the 
radial AdS coordinate gets too close to the boundary. If z is the radial AdS^ coordinate 
we say that z > ec- Here ec is a short distance UV cutoff from the Wilson loop perspective 
or an IR scale from the amplitude point of view. Thus we are cutting off the z integral 
in a way that depends on the solution. However, this dependence on the solution can be 
expressed in terms of the distances between the cusps (it depends on distances, not only 
on cross ratios). We explain this in detail in the appendix. We find that 

A = Adiv + ^BDS-like + ^periods + ^free + COUStaut (4.11) 

where Adiv contains the ec dependent divergent terms. AeDS-iike is finite but it depends 
on physical distances and not just the cross ratios. Its explicit form is given in appendix 
A. It obeys the anomalous conformal symmetry Ward identities derived in [|3^. It has 
become conventional to define a remainder function by subtracting the specific function 



Abds given in [35|, which has the functional form of the one loop result. Thus one defines 



a remainder function via 

R = ~ {A — Adiv — Abds) = —{ABDS-iike — Abds) — ^periods — ^free + constant 
R = ^ ~^ ^-^^2(1 - Ui)j - \Z\'^ - Afree + constant 

(4.12) 

In the large \Z\ limit, which is the coUinear limit, one can check using (|3.9| ) that the 
first term in R, involving Ui, cancels the |Zp term, up to a constant. Thus the remainder 
function goes to a constant in the large \Z\ limit, since A free — )■ in this limit. As argued 
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35(1, the Abds term gives the expected result in the coUinear hmit. Thus, the remainder 



function should be finite. This is a simple check of our results. 

Let us summarize the final result. The objective is to find the remainder function 
i? as a function of the spacetime cross ratios Ui. For this purpose, we compute fx from 
the Ui via (|3.14| ) and we introduce an auxiliary complex parameter Z. We then solve the 
integral equations (|3.6|) and determine the functions e(6'), e{9). These functions depend on 
the parameters \Z\ and fi. We then solve for the spacetime cross ratios Ui (or equivalently 
bi ( p.lO|) ) using (|3.7|), (p78|). In principle, one would be able to invert that relation and 
compute Z and ^ in terms of the spacetime cross ratios Ui. Finally, one can compute A free 
as a function of |Z| and fj, using ( [4.10| ) and insert the result in ( [4.12| ). In practice one can 
numerically compute the functions e, e, the Ui and the area as a function of |Z| and 



5. Analitic solution of the integral equation 

In the following we will compute the area from the integral equation in the par- 
ticular limit [/ — 7- 0. This corresponds to the special kinematical configurations with 
ui = U2 = Us = u. From the point of view of the TBA equations this limit corresponds to 
a high temper ature/conformal limit and has been analyzed for instance in |]3^,0 (see also 
1 38| , ^ , |40[| ) to which we refer the reader for the details. We start by recalling the general 
result in these references and then we focus on the specific case at hand. Consider the 
standard form of TBA equations 

/dt' 
—Kabit - t') log(l + AbC-^''^*')) (5.1) 

where the indices a, h run over different species, nia denote their masses and \a their 
chemical potential. Then, the free energy of the system can be written as 

/"(^tcoshtlog(l + Aae-^"W) (5.2) 
^ 2ti J 



In turns our that the free energy can be exactly computed in this limit [p6| , |37| , p8ip9 



F{ma ^ 0) = -^Y^CxM. Cx{x) = \ (\o^x\o^{\ + ^) - 2L.2 
The parameters Xa are the solution of the system of equations 



(5.3) 



Xa 

b 



n(l + - , Nab = -— dtKabit) (5.4) 



Xf, I Zn 



oo 



Note that the final result does not depend much on the details of the Kernels (only their 
integrals) 
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5.1. Particular case at hand 

The integral equations ( |3.6|) are exactly in the standard form ' that we can 

apply the above results straightforwardly. We have three species of particles, that we call 
ei, 62 and e (with ei = £2 = e). The mass of each excitation is i 

mg = V2\Z\, rrie^ = rue^ = \Z\ (5.5) 



The chemical potential for each species is Ag = 1, A^^ = n and Aej = From the 

particular kernels in our equation we can easily compute Nab- We find 



iVgg = 1, iVg,, = iV,^g = 1, iV,,,^. = 1 (5.6) 
Pluggin this in (|5.4|), we can solve for xg and — x^^ = x^^, we find 

= xg = 1 + ^2/3 + ^-2/3^ e' = x, = ^1/3 + ^-1/3 (5.7) 

The free energy is simply 

Using the explicit expression for C and identities involving di-logarithms one can show the 
very simple result 

-F = A/,ee = ^-^</>', ^i = e'^ (5.9) 

Several comments are in order. First, is it possible to solve numerically the integral equa- 
tion for small values of Z (for instance Z = 0.01). After several iterations we find numerical 
results that are in perfect agreement with the analytic solution ( [5.9| ) , see for instance the 
following figure 



^ Quite interestingly, note that the mass of these three excitations agree (up to the propor- 
tionality factor \Z\) with the masses of the three fluctuations around the classical string solution 



describing a single cusp (or four cusps) in AdSs, see [41| . 
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A 




Fig. 14: Analytic (continuous line) vs. numeric (points) solutions of the integral 
equation for the case Z — 0.01 and /x = e*"^, with running from to tt. 



In particular, this shows that the numerical solutions can be quite reliable. 

Second, we see that the values of the regularized area/free energy for = ±1 is simply 
given by ± ^ . This values of correspond to u = 1/4 and u = 1 respectively, and the value 
for the free energy is in very good agreement with what we expect from the numerical 
analysis of appendix B. Recall that /i = ±1 corresponds to the two kinds of regular 
hexagons. 

Third, the free energy can be expressed in terms of u by using 

1 1 — 3u , , 

fi + ^ =2cos(p= 5.10 

where (p is not necessarily real. We are interested in the whole region u > 0. We can cover 
this region by going along the following contour in the //—plane 



Imn 






u=l / 


^u=l/4 


u=0 




U=oo 


Rep, 



Fig. 15: Different values of u as we follow the path shown by the figure in the 
/i— plane. For u > 1/4 fx is a phase. 



The substraction of 2Apentagon in appendix B is needed since the free energy was defined in 
such a way that is zero for very large values of \Z\. 
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As n is real and very large, u is very close to zero. As we approach n = 1, u grows 
until u = 1/4 aX fx = 1. For u > 1/4, /i is a phase, with /i = — 1 at w = 1 and u becoming 
very large as fx approaches —i. 

The final answer for the remainder function in this regime is then obtained by adding 
up all the contributions 

1 . 5? „ 1 

(5.11) 



TT 1,2 3 



R{u, u,u) = -^ + ^(f' + - (log^ u + 2Lz2(l - u)) 
6 Stt 8 ' 



u 



4cos2(0/3) 

This is the remainder function for the scattering of six gluons at strong coupling in the 
particular kinematical configuration in which all the cross-ratios ui coincide. 




log(nl 
log(lO) 



Fig. 16: Remainder function at strong coupling for 



Ul = U2 = U2 



5.2. A curious observation 



The remainder function at two loops for the case at hand was extensively considered in 
|4^ , based on previous work [13:1^,43|. In that paper the remainder function was computed 
numerically for several values of u (in particular, u = 1/9, 1/4, 1,3.83 and u = 100). We 
could try to fit their numerics by a function with some arbitrary coefficients but the general 



structure of ( 5.11 ) , similar to what it was done in [44] for the case of the octagon in AdSs- 
More precisely, we consider a function of the form 

Rc„c,,cs = Cl + ^'/''^ + C2^ (log^ U + 2Li2{l - U)) + C3 (5.12) 

The idea is that the high temperature of the TBA equations are not too sensitive to the 
precise form of the kernels, so that perhaps this functional form holds for all values of the 
coupling. In addition, in our computation, the two terms arise in a somewhat independent 
fashion, so we have given us the freedom to change their relative coefficient^!. It turns 
out that for certain values of c's, namely ci ~ 12.2, C2 ~ 11.4, C3 ~ —9.1, we get a quite 
good approximation of the two loops result, see the following figure 



Note that we obtain that Ci/c2 ~ 1.07, which is close to one, so perhaps we should not change 
the relative coefficients of the two terms! 
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Fig. 17: -Rci,c2,c3 (continuous line) vs. numerical values (points) for the remainder 
function at two loops 



However, the asymptotic value for as w — ?■ oo seems slightly higher than the value 
computed in |^ . It would be interesting to try to understand better the two loops result in 
this particular regime. In particular, it would be interesting to find its analytical expression 
and compare it with the strong coupling result. 

Comparing the two loops (figure 8 of ) and strong coupling results , we see that 
many qualitative features are very similar and are probably independent on the coupling. 
For instance, the apparently universal behavior R ^ h{X) log^ u for small u may be simple 
to understand and h{X) may be computable for all A. 

After the first version of this paper was published, an analytic expression for the two 
loop result was found [^. The function looks very similar to the strong coupling result, 
but it cannot be fitted by an ansatz of the form (|5.12|). 



6. Conclusions 

In this paper we have studied minimal surfaces in AdS^. These are the classical 
solutions for strings moving in AdS^. We have shown that the problem can be mapped 
into an SU{4) Hitchin problem. More precisely, it is a Z4 projection of an SU{2, 2) Hitchin 
problem. An important observation is the existence of a spin four holomorphic current 
P{z) on the worldsheet. For surfaces that end on a polygonal contour, P is a polynomial 
whose degree depends on the number of sides of the polygon. In the Hitchin language 
we have an essential singularity at ^ = 00. This singularity has its associated Stokes 
phenomenon. This Stokes phenomenon is important for the emergence of the polygon on 
the boundary of AdS. In fact, the shape of the polygon is related to the Stokes matrices 
of the solutions of the fiat Hitchin connection. The shape of the polygon is encoded in the 
coefficients of the polynomial as certain phases which correspond to other moduli of the 
Hitchin system. 
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The integrability of the Hitchin equations can be used to analyze this problem. We 
introduced the parameter ( and studied the problem as a function of (. When ^ is a phase, 
we have the ordinary physical problem, but with complex C it is some sort of deformation 
of the physical problem. The nice feature is that for C — or — )■ oo one can solve it 
using a simple WKB approximation. This solution displays a Stokes-like phenomenon but 
now in the ( plane. The Stokes "factors" can be found explicitly by analyzing the small 
C problem. Given that we have an analytic function of ( (analytic away for ( =, oo) with 
given Stokes-like behavior at zero or infinity, we can recover the function. This is done 
by first formulating a Riemann-Hilbert problem and then writing its associated integral 
equation. 

Surprisingly, this integral equation ( |3.6| ) has the form of a TBA equation with rather 
simple kernels. In fact, they are the TBA equations for an A3 (or Z4) theory discussed in 
Il30| , p8| . The area coincides with the free energy of the TBA system. 

The TBA equation can be solved rather simply by a numerical iterative process. The 
TBA equation depends on two parameters \Z\ and fj, (|3.6|) . The non-trivial part of the 



area, which is the free energy of the TBA, also depends on \Z\ and n ( [4.101 ). Once the 
solution is found we can determine the spacetime cross ratios, wi, M2, us, in terms of Z 
and fx (|3.7| ) ( P^) . This leads to the final result for the remainder function (|4.12|) . 

There is a special set of cross ratios ui = U2 = for which the problem simplifies. 
This corresponds to the high temperature limit of the TBA. This limit was studied in 
1^,^. Using these results we could find an explicit, and rather simple, analytic expression 
for the full answer (|5.1ip . 

For us, the connection to the TBA system was an "experimental" observation. It 
would be nice to find a simple argument leading more directly to this TBA equation. It 
is interesting that the TBA that we found corresponds to the A3 theory ||3^,38|, which is 
associated to S'f/(4), the group associated to the Hitchin equation^. 

Extending the techniques used in this paper to the case of more gluons should be 
straightforward. We expect that we will end up with a more complicated pattern of 
discontinuity lines in the ( plane. We expect that the lines will organize in groups of four 
lines at right angles, but with different groups rotated with respect to each other. The 



Note that the spectrum of particles in the A3 integrable theory |3C,3S], two particles of 
mass m and one particle of mass \/2m coincides precisely with the spectrum of particles for small 



fluctuations around the string worldsheet describing a single cusp, analyzed in [41|. 
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TBA will have a similar structure, probably with one unknown function per group of four 
lines. But the final area is still given by the free energy of this TBA system. 

One would also like to generalize these results to the full quantum problem. It is 
natural to imagine that perhaps the full quantum answer is also given by a TBA equation 
involving the same particles but with a different kernel which includes all quantum correc- 
tions. Or perhaps we have a similar kernel but a different mapping between the parameters 
Z and \x of the integral equations with the spacetime cross ratios. A natural next step 
would be to extend the Pohlmeyer reduction and the relation to the Hitchin equations for 
the full PS {2, 2|4) supercoset theory. This should make it clear how to compute N'^MHV 
amplitudes. 

The large temperature limit of the TBA calculation is insensitive to the precise form 
of the kernel. What is important is the total integral of the kernel (|5.4| ). Thus, one can 
speculate that perhaps for all values of the coupling one has a form similar to the one 
we found, but perhaps with a different coefficient ( |5.12| ). In fact, if one allows for general 
coefficients, the functional form that we found is very close. 
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Appendix A. Regular izat ion of the area 

In this appendix we regularize the area using a physical cutoff which is appropriate 
to regularize the Wilson loop or the amplitude. The idea is very similar to the one used 
in For the sake of completeness we give a brief description. The story is relatively 
simple in the case that the number of gluons is not a multiple of four, n ^ 4k. We will 
treat only this case here, since we mainly deal with the n — 6 case which certainly obeys 
this condition 111. 



The case where n = Ak was treated explicitly in [p5| for AdSz kinematics. We expect similar 
results for generic AdS^ kinematics. 
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The area can be written as 



-A — Afj-ee ^periods "l~ ^cutoff 

^periods =2 / d^z{PPy/^-2 I d^w = 2 I I d^w 



(A.l) 



^cutoff =2 / d w 

Where A free is the non-trivial function that is computed as the free energy of the TBA 
equation (|4.9|) for the case of n = 6. In the definition of Aperiods we have used that the 
polynomial P defines a surface E which is the w plane with the appropriate cuts. We 
have introduced Eq as a reference surface with a single branch point at the origin and the 
same structure at infinity. In order to subtract the two it is very important to regularize 
both computations with a cutoff in the w plane, such as = A ^ 1, but not a cutoff 
in the z plane. Using this cutoff for defining the difference in Aperiods we find that it 
can be explicitly evaluated in terms of periods on the Riemann surface x'^ = P{z). It 
has the form of Ze^Zrm — Z^.Ze^ where e^, rrii is a basis of electric and magnetic cycles 
(i.e. a basis of cycles which are symplectically normalized). For the particular case of the 
hexagon, Aperiods = , where Z is defined in (|3.1| ), see (^J]). Acutoff is imposing the 
physical cutoff which requires that the AdS radial coordinate should be larger that 
l/cc- Here ec is a distance scale providing a UV (IR) cutoff in the Wilson loop (amplitude) 
interpretation. 

In order to evaluate Acutoff we can take A = and replace the w surface by a simple 
surface with all the branch points at the origin. Our task is now to evaluate AcutofF- In 
order to evaluate it we will need to use a couple of properties of the solutions discussed in 
section 2.4. Near each cusp we have a behavior which goes roughly like 

A „.A ( 2Re\w\ „^ 2Im\w\ -2Re\w'\ _ -21 



y^^^yf X (^e'^^^'^J, or e'^'^^'^J, or e'^^^^'^J, or e-^'"^^"'Jj (A.2) 

with a particular choice at each cusp (we have rotated the w plane so simplify these 
expressions). The position of the cusp on the boundary is determined by the null vec- 
tor yi {yi-Ui = 0). Recall that we had a basis of small solutions s^, normalized so that 
det{si, Si_|_2, ■Si+s) = 1- The behavior of the spacetime solution near the cusp i could 
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be found as = qfsfV^ Si-\-i. Which, together with the normahzation condition on the 
imphes 



yi.yi+2 = 1 



(A.3) 



Let us say that at cusp i we have Yi ~ y^e 



,e[w]^ Then the spacetime cutoff has the form 




+ 2Re[w] 



(A.4) 



which gives a hne in the w plane, 2Re[w] = — log — log y^ . We have a similar relation on 
the cusp Y-'^2j where the line is now at 2Re[w] = — (— log[ec] — logyj^2)- We first imagine 
some reference lines at 2Re[w] = ilogCc which are then displaced by small distances 
proportional to 5i = —logy^. It is convenient to relate the 5i to the spacetime distances. 
This is done as follows. We know that yi.yij^2 = 1- At the same time we know that the 
spacetime distances can be expressed as (i? , , ~ This implies that 



where we used (|A.3| ). This equation, together with Si+n = can be used to solve for all 
the 5i in terms of the logs of the nearby distances li. Here it is important that n is not a 
multiple of four. 




5i + 5,+2 = ii = log 



'i,i+2 



(A.5) 




branch 



cut 




Fig. 18: Here we show the computation of Acutoff- The w surface So has a 
single branch point at the origin. We show only a portion of this surface. Here 
L = — log Ec and 5i = — log . The total area is obviously the sum of the areas of 
the various rectangles. Note that L ^ 5i as we take the cutoff away. 
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We can then compute the area by summing the areas of various rectangles in fig. 18. 
We get 

1 



^cutoff 
^cutoff 



2 



+ ^BDS-like 

^div = - 2^ ( - log + ^ 



(A.6) 



^BDS-like 

We can rewrite this as 



^BDS- 



like 



n / 2K 

i=l \ k=0 



^fc+l 




l+2k 



(-1) 



fc+1 



for n = 4K + 2 



for n = + 1 , 4K + 3 



(A.7) 

Here Adiv is the standard divergent term that is present for hght-hke Wilson loops or 
for amplitudes [|6.5^,|8|. On the other hand AeDS-iike is a finite term which obeys the 



anomalous Ward identities of broken conformal invariance ||3^ . An easy way to check this 



is to note that if we shift each log of a distance logd^^ — )■ logd'f^ + + e, , then 



dAsDS- 



ike 



dA 



BDS 



dti 



|e=0 



dei 



|e=0 



1 d'^ 

^ "■i,i+2"i,i-2 



where Abds is the expression written in which is the one loop result in the weak 
coupling expansion, except that here it comes with a different overall factor, due to the 
different value of the cusp anomalous dimension. When we act with the dual conformal 
generators we see that we need to add with these derivatives. Thus if we have the same 
derivatives as the BDS answer, then it obeys the same equations as the BDS answer. Note 
that if n is not a multiple of four we cannot make any cross ratios from the distances (ii,i+2- 
This implies that there is a unique solution of the anomalous ward identity which can be 
written only in terms of (ii,i+2- This unique solution is ABDS-iike- 

In order to find the difference between Abbs ~ ^BDS-iike we simply start with any 
non- nearby distance, dij, appearing in Abds and we write it in terms of the unique cross 
ratio Cij which involves dij and nearby distances (ii,i+2- 
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Thus, if one is interested in the remainder function, which is conventionaUy defined 
in terms of what we have to add to Abbs to get the fuU answer, then it can be expressed 
as foUows 

~A = - Adiv - Abds + R 

(A.9) 

R = - (ABDS-like - ^BDS) - ^periods - ^free + COUStaut 

where the constant is independent of the kinematics. 

In the AdSs hmit these formulas become the expressions for n/2 odd in (Formula 

5.8). 

In this paper we are mostly interested in the case n = 6. In this case, the BDS ansatz 
readJi^ 
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^BDS = ^ 



1=1 



1^2 ^2 1 ^2 1 / ^2 /^2 

_ lo ^'"^"^^ lo — lo ^ '^2,2+3 _|_ jL I ^ _ ■^i,i+2-^z+3,i+ 5 

A ^2 ™2 1 ™2 Q I ™2 ™2 

^ -^i.i+S "^i.i+S -^2+1,1+4 ° \ -^1,1+3-^ i+2,i+5 



(A.IO) 

We find that 

Abds - ^BDS-like = ^ (^log^ Ui + ^Lz2(l - Ui)j (A. 11) 

Hence, the remainder function is obtained by adding the above difference to the regularized 
discussed in the body of the paper. 



Appendix B. A new class of regular solutions 

In the following we will consider a generalization of the regular polygon solutions 
considered in [^5|] . These are simple solutions for which the area can be evaluated directly. 



Thus, they provide special cases which can be compared to the general results we obtained 
above. We consider solutions that possesses Zn symmetry. In other words, we have 



V{ze'-,ze-'-) = R%Y''{z,z) (B.l) 

where R a rotation in space-time, R^ = 1. These solutions are associated to regular 
polygons. Using the expression (|2.1C1| ) for the holomorphic function P{z), this implies 



; 2tt . 



P{ze'^)=e-^'^P{z) (B.2) 



Note that what we are calling Abds is minus the expression in [^5| (up to a factor) because 
of the minus sign in ( |2.8|) . 
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The simplest holomorphic function satisfying this condition is a monic polynomial P{z) = 
Pqj^ i^g^gg ^ ^ ^ "^^g h ave n even. In the general case of AdS^, the simplest 
elementary solution would correspond to P{z) = z, which would correspond to a regular 
pentagon. 

As seen in it is convenient to use (2,2) signature to construct the pentagon 



because we can then choose all non-consecutive cusps to be spacelike separated. Further- 
more, we will see that it can be actually embedded into ^^5*4, very much like the regular 



polygons studied in could be embedded in AdSs. 

We start by considering the embedding coordinates of AdS^ in the case of (2, 2) 
signature. 

-y^2 - Y-i - Yq + Yi + Yi + Yi = -1 (B.3) 

the restriction to AdS^ corresponds to setting Y3 = 0. Note that this is an AdS4^ space 
with two times, whose boundary has two time and one space directions. We can instead 
use ordinary AdS^ Poincare coordinates, which correspond to ^ = Y-2 + Ys = Y'_2, 
(^1,^2) = ^'^y^'J'"'* and (xi,X2) = ^^y^^^ ■ '^^^ regular polygons to be considered in this 
appendix can be embedded into this AdS4^, hence, they correspond to contours in the 
(xi, X2) plane and in the (ti,t2) plane, such that 

xj + xj- 4-4 = 1 (B.4) 

This is a surface inside i?^'^ with two time and one space directions. This surface is the 
boundary of an AdS^ subspace of AdS^ where the surface is residing. In addition, we 
require each segment to be light-like, so the lengths in the a;— plane and t— plane coincide. 
Finally, we consider only polygons such that at each cusp the momentum transfer is space- 
like. This can be achieved by requiring that the angle between consecutive segments in the 
x— plane is larger than the angle in the t— plane. A prototypical example of such contour 
is the regular pentagon 

(t,,t2) 





Fig. 19: Regular pentagon in (2, 2) signature. The lengths of each segment in the 
x— plane and t— plane coincide and the angles are always smaller in the t— plane. 
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B.l. Strings on AdS^i 

Next, let us briefly consider the reduction of strings on AdS^. In this case Bi and B2 
are non trivial and we have vi and f 2 correspondingly. From now on, let us focus on the 
case of (3, 1) signature and parametrize 

vi = iP{zy/^ cos (/3/2), V2 = P{zy/^ sin(/3/2) (B.5) 

where P is real. Later on, we will then continue to (2, 2) signature by simply taking /3 — > i^. 
We would like to express the physical field /3 purely in terms of space-time quantities. For 
that, let us recall the following expressions 

d^Y = dadY + v^Bj, d^Y = BadY + v^Bi (B.6) 

Viewing this as an expression for ^ we can compute 

-v^v^ = \P[z)\ cos/3 = e'^dada - d^Y.B'^Y 
In addition, with a little bit of effort, one can show 

e°'\P{z)\ sin /3 = eabcdeY^dY^dY^d^Y'^d^Y^ (B.8) 

From this relation we can immediately see that (3 is real and also that (3 ^ ij3 is the correct 
prescription in order to continue (3, 1) signature into (2, 2) signature. 

B.2. Boundary conditions at zero 

Let us now focus on the case of (2, 2) signature. In that case, the equations of motion 
and Virasoro constraints imply the following equations for d = a — | log \ P{z)\ and ^ 

dda = e^- e'^ cosh /3, dd^ = e"" sinh ^ (B.9) 

where we have written these equations in the w— plane, defined by dw = P{z)^/^dz.A.s it 
can be seen, a and (3 are massive fields that decay exponentially at infinity. 

In order to understand the boundary conditions as we approach one of the zeros of 
P{z), we use polar coordinates z = pe^'^, z = e~'^'^ and consider the expansion of a general 
solution around p = 

Y_2^1 + J2^i~'\v)p', = i = -1,0, 1,2. (B.IO) 
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(B.7) 



By considering the equations of motion order by order in p, the functions ^^^\(p) are 
determined, up to integration constants at each order. In addition, the requirement = 
— 1 and the Virasoro constraints impose some conditions on these integration constants. 
When considering the case of the regular pentagon, we have an additional constraint, 
coming from the Z5 symmetry of the contours, see fig. 19 . More precisely, we can introduce 
complex spacetime coordinates X = Yi + iY2, T = Y-i + iYq and we then require 

X{ze'^,ze-'^) = e'^''-X{z,z) , T{ze'^ ,ze-'^) = e'^'''T{z,z) (B.ll) 

Where r^; and are two integers which parametrize the spacetime Z„ rotation matrix. 
For the pentagon case, n = = 1 and = 2, see fig. 19. Once all the constraints are 
worked out, we can compute the expansions of a and /3 around the origin. As expected, a 
is regular, while 

,^ e^dada-d'^Y.d'^Y 1 

Hence (3 has a logarithmic divergence as we approach p — )■ 0. A similar argument for (3, 1) 
signature, shows that we should expect some winding in (3 as we approach the origin. Since 
cosh /? > 1 and cos /3 < 1 it is reasonable to expect such behaviors for /3 from ( |B.12| ) 111. 

Hence we expect the regular pentagon in (2, 2) signature to correspond to two fields 
a and (3 with logarithmic divergences at the origin and exponential decay at infinity. 
Furthermore, it is reasonable to assume that such solutions are rotationally symmetric. 
Namely a and (3 are functions of |w| = p only with the following boundary conditions. 

2 4 

= -3 log p + c« + /3(p) = --log p + C/3 + p^O 

5 5 (B.13) 

d(p) ^ Ko(2y2p), /3(p) ^Ko(2p), p^oo 



15 



More precisely, the left hand side of ( [B.12 ) will have a general expansion of the form 



r , „ „, e°'dada — d^Y.d^Y c\ + C2Z + + ... 
[cosh/3 or cos/3] = \Q2Y,d'^Y\ " ~p 

For (2,2) signature, as cosh/3 > 1, we expect ci 7^ 0. On the other hand, for (2,2) signature, as 
cos/3 < 1 we expect ci = and C2Z + C2Z = 003(99 + 7), resulting in non-trivial winding for /3. 
The winding for /3 in (3, 1) signature is consistent with the assumption of /^Dorn as defined by [23|, 



being regular everywhere. In | ]23[ | P is factorized in two polynomials P{z) = V-{-{z)v-{z) and the 
relation between our /3 and the one used in [23| is /?ours = /3Dorn — f log ^ • 
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The constants Cq, and cp have to be chosen in such a way that the solutions decay at infinity. 
We find numericaUy that such solutions indeed exist and their exponential decaying is in 
perfect agreement with the expectations. Solving numerically the shooting problem, we 
find Ca ~ —0.28, C/3 ~ —0.38. These values are consistent with a crude version of the 
analysis of 0. Once the numerical solution is found, we can compute its regularized 
area. We find the approximate result 

^pentagon = 2 ^(e" - l)d^W ^ 1.18 (B.14) 

with an estimated error of a few percent. Note that /3 — t- — /? is a symmetry of the equations 
and we get another solution where (3 log p at the origin. These two pentagons differ 

by a spacetime parity operation^. For a more general polynomial we can choose either 
boundary condition, — )■ ±| logp, at each zero of P{z). 

B.3. Other regular solutions 

Now we can imagine that P{z) has two zeros, in each of which (3 has boundary 
conditions (5^. As we bring the two zeros together there are two distinct situations: 
1.- We bring together two zeros where the boundary conditions are (3^ and In this 
limit, we expect to recover the usual regular hexagon, which can be embedded into AdS^ 
and for which /3 = 0. In terms of contours in (2, 2) signature, its boundary looks like 

(X,,X2) (ti,t2) 




2,4,6 



1,3,5 



Fig. 20: Usual regular hexagon in (2, 2) signature. As it can be embedded into 
AdSs,, since one of the time coordinates vanishes identically. The cross ratios are 
U\ = U2 = = 1. 



However, in |5(][ | it was assumed that /3 is regular at p — > 0, while we assume a logarithmic 
divergence (|B.13|)r. 
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The fact that spacetime parity changes the sign of /3 can be seen from 
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The regularized area corresponding to this hexagon was computed in [25|, with the 
result 

2.- We bring together two zeros where the boundary conditions are both /?+ (or both (3~). 
In this limit, we expect to recover a different regular hexagon, still with symmetry, but 
which cannot be embedded in AdSs because (3^0. Furthermore, the field (3 corresponding 
to this solution has a logarithmic singularity at the origin. In terms of contours in (2, 2) 
signature, its boundary looks like 



(X,,X2) 



(tl,t2) 






Fig. 21: Alternative regular hexagon in (2,2) signature. It still posses a Zq 
symmetry but cannot be embedded into AdSs. The cross ratios are u\ = U2 = 

U3 



1 

4 • 



In the w— plane we expect the following boundary conditions at the origin 

2 4 
a{p) = -- log p + Co, + ..., (3{p) = --\ogp + cp + ..., p-^0 (B.16) 

The condition for (3 is deduced by looking at ( |B.71 ) near the origin and demanding that 
the right hand side is approximately constant near the origin. Solving numerically the 
shooting problem and computing the approximate area, we find 

<tagon ^ 2.85 (B.17) 

The value of the cross-ratios is = 1 for the first, usual, regular hexagon and Ui = 1/4 
for the second kind of regular hexagon. It is difficult to compare these individual areas 
to the result of the integral equation due to possible overall additive constants. However 
we can compute the hexagon areas minus twice the pentagon area, which corresponds to 
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the difference in the regularized area when the two zeros are brought together from a very 
large distance. From our numerical estimates we find 

TT TT 

^hexagon 2Apentagon ~ -0.52 ^ --, ^hetagon ~ ^Apentagon ~ 0.50 - (B.18) 

The values ±7r/6 are the exact results obtained using the integral equation. 

For a generic number of zeros of P{z) we expect a similar situation. Given a polygon 
of n-sides, we can draw a regular polygon such that we need to perform a rotation by 
u = e^2^p/^ in order to go from one cusp to the next one, with p = 1, 2, .... Hence, we will 
have a family of inequivalent polygons labeled by two integers {pxiPt)-, drawing a regular 
polygon with p^^ in the x— plane and with pt in the t— plane, see ( [B.ll|) . For instance, in 
the following figure we see different choices p = 1, 2, 3, 4 for n = 9. 



(■) (>] (c) (a) 




Fig. 22: Shape of different regular nonagons, with p = 1, ...,4 from left to right. 
The full polygon in (2,2) signature is given by choosing a pair of these with px < pt- 



Furthermore, the requirement of space-like momentum transfer can be achieved by 
imposing Px < Pt ^ [f]- Finally, let us mention that such polygons can be embedded 
into AdS^ and the distance between adjacent cusps can be made light-like. Centering the 
above polygons at the origin we can rescale their relative size so that the segments in the 
a;-plane and the t-plane have the same length. Then the pair of regular polygons produces 
a null-polygon in R^'^. 

It would be interesting to understand these solutions better. It is likely that they cor- 
respond to radially symmetric solutions. One could attempt to find the precise boundary 
conditions for /3 in each case. 

Appendix C. On the appearance of poles in the TBA equation 

Although in the initial large \Z\ region the Xi[Q had no zeroes or poles in C, they 
may appear at smaller \Z\. This is a familiar phenomenon in TBA-like equations p3| , p2| . 
Poles and zeroes appear in a very specific manner. Lets' start from a situation without 
poles or zeroes, and for simplicity consider some function with a discontinuity = 
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+ y) across a cut £. In general, both the X'^, X~ which define X on the two sides 
will be nice analytic functions across the cut. A'+ is supposed to have no zeroes or poles 
above the cut, but it may have either below the cut. It ie easy to see it has to be a zero: 
X~ and y are not supposed to have zeroes or poles there, all can happen is that 1 + 3^ is 
zero at some point Co- 

As we vary the parameters in the integral equation, the point (o defined by 3^ [Co] = — 1 
may move across the cut. In that case, we do not expect the zero of X'^ to suddenly 
disappear, and a pole of X~ to suddenly appear if we follow our solution continuously. 
That would happen, though, if we did not modify the integral equation. The reason is 
clear: the integral equation is of the form logX — ■ ■ ■ + Ki * log(l + y): the log has a 
singularity at Co, which we are passing across the cut. 

The correct way to continuously follow a solution is to accept that X~^ should have 
a zero at Co above the cut. The integral equation is easily modified to account for it, 
schematically as 

\ogX=--. + log(C - Co) + Ki * log(l + y) 

The extra term actually coincides with the difference between the convolution Ke*log{l+y) 
along paths just above and just below Co- Notice that the location of the zero is fixed by 
3^ [Co] = —1, in order for X~ not to have the zero there. 

The expression for the TBA free energy is similarly modified by the addition of a 
schematic term Z/(o + ZQq. In TBA these extra terms have the interpretation of an 
excitation added to the ground state. In our specific setup, the location of possible zeroes 
and poles is determined by equations of the form ttjCo] = oo from the cuts and 
equations of the form 6j[Co] = "A* or 6i[Co] = —fJ'~^- It is easy to see that if, say, bi = —ji 
then either = —jjT^ or = —jjT^^ and one of the cross-ratios equals 1. The 
converse is also true. It seems possible that as we follow a path towards (3, 1) signature 
such phenomena will occur. 
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